Conservative upper limits on WIMP annihilation cross section from Fermi-LAT 7-rays 
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The spectrum of an isotropic extragalactic 7-ray background (EGB) has been measured by the 
Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subtrac- 
tion of unresolved point sources and extragalactic diffuse processes, which could explain from 30% 
to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due 
to the annihilation of dark matter (DM) particles in the Galactic halo, we obtain stringent upper 
limits on their annihilation cross section. Severe bounds on a possible Sommerfeld enhancement 
of the annihilation cross section are set as well. Finally, we consider models for DM annihilation 
depending on the inverse of the velocity and associate the EGBs to photons arising from the anni- 
hilation of DM in primordial halos. Given our choices for the EGB and the minimal DM modelling, 
the derived upper bounds are claimed to be conservative. 
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I. INTRODUCTION 



The indirect search for dark matter (DM) through its 
annihilation products in rare charged cosmic rays (CRs) 
and in multi-wavelength channels requires very accurate 
measurements and an unambiguous estimation of all the 
possible backgrounds to the DM signal. In the last years, 
dedicated experiments have provided unprecedented re- 
sults by extending the energy ranges of the measured 
cosmic species as well as the precision of the data 

Further data are expected by the Fermi-LAT and 
Pamela on-going missions, and by the AMS-02 experi- 
ment on board the International Space Station. From 
the theoretical side, many efforts have been addressed to 
a better and increasingly detailed modellization of the 
astrophysical processes which shape, at different levels, 
the observed fluxes. Data from cosmic antiprotons Q 
have been shown to be compatible with the standard pro- 
duction from CRs impinging on the interstellar gas [1]. 
The anomalous increasing positron fraction measured by 
Pamela [U, H and confirmed by Fermi-LAT may be 
explained by emission from near pulsars over-imposed 
to a standard CR population [l^, [ll|. Alternatively, a 
DM component with very high cross section or sources 
concentration has been invoked Unprecedented 
7-ray measurements by Fermi-LAT have boosted inter- 
pretation of diffused and point sources emission in terms 
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of exotic components from DM annihilation in the halo of 
the Milky Way, in ext rag alactic near objects or in cosmo- 
logical structures 0, [la - iTsI . The very signature would 
be the monochromatic line, which nevertheless provides 
tiny signal on a remarkable background [ISy. 
The high latitude 7-ray emission measured by Fermi-LAT 
[t'], given its reduced contamination by galactic sources, 
can be a powerful tool to set limits on the contribution 
of DM to the measured flux. The data are indeed the 
result of a non trivial subtraction procedure and show a 
high isotropic feature. 

The aim of the present research is to set conservative 
upper limits on the galactic weakly interacting massive 
particle (WIMP) DM annihilation cross section into 7- 
rays. Several upper limits have been obtained through 
different and complementary indirect research means 
d [H HO? -[13 ■ However, it is usually not straight- 
forward to compare these results, given the model de- 
pendence, the different assumptions on the astrophysical 
backgrounds, and the theoretical uncertainties. We will 
confront the 7-rays coming both from the DM halo and 
high-redshift protohalos with the background observed 
by Fermi-LAT at high latitudes. The conservative ap- 
proach is achieved - in addition to prudent assumptions 
on the particle physics model and DM distribution in the 
Galaxy - through the comparison of the putative DM sig- 
nal with a high latitude diffuse emission spectrum (i.e. 
EGB) obtained with minimal subtractions of known un- 
resolved sources. 

Our paper proceeds as follows. In Sect. 2 we discuss the 
possible contributions to the high latitude 7-ray emission 
from unresolved point sources and truly diffuse processes. 
We subtract the non-negligible fluxes to the Fermi-LAT 
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data and draw two possible scenarios for the high lati- 
tude emission. In Sect. 3 we derive conservative upper 
limits to the DM annihilation cross section by identifying 
the residual 7-ray flux with 7-rays from DM annihilation 
in the galactic halo and in primordial DM small halos 
at high redshift. In the latter case, we study models in 
which the DM annihilation cross section has an explicit 
dependence on the inverse of the velocity. We discuss 
also a possible Sommerfeld enhancement of the annihila- 
tion cross section and derive limits on its amplitude. In 
Sect. 4 we draw our conclusions. 



II. THE EXTRAGALACTIC 7-RAY 
BACKGROUND 



A diffuse 7-ray emission has been measured by the 
Fermi-LAT detector at high latitudes {\b\ > 10°) The 
spectrum has been obtained after the subtraction from 
the data of the sources resolved by the telescope, the (in- 
deed model dependent) diffuse galactic emission, the CR 
background in the detector and the solar 7-ray emission. 
The resulting flux decreases with a power law of the pho- 
ton energy with spectral index 2.41 ± 0.05. It shows a 
highly isotropic sky distribution and is generically classi- 
fied as an extragalactic 7-ray background (EGB). 
The 1451 sources listed in the First Fermi-LAT catalog 
(IFGL) 28] represent the best-resolved survey of the sky 
in the 100 MeV to 100 GeV energy range. For each low- 
flux source there may be a large number of unresolved 
point sources which have not been detected because of 
selection effects, or too low emission. Most of the unas- 
sociated high latitude sources are blazars, a class of Ac- 
tive Galactic Nuclei (AGNs), and their pile to the EGB 
with the largest flux 29]. Galactic resolved pulsars and 
Milli-Second Pulsars (MSPs) represent the second largest 
population in the Fermi-LAT catalog [H, [s^l and they 
are expected to contribute significantly to the putative 
EGB. A non- negligible 7-ray flux seems to be guaranteed 
by unresolved normal star- forming galaxies [31|. Ultra- 
high energy CRs (UHECRs) may induce secondary elec- 
tromagnetic cascades, originating neutrinos and 7-rays 
at Fermi-LAT energies j32i ]. Contributions from unre- 
solved blazars and MSPs are believed to contribute at 
least few percent to the Fermi-LAT EGB, while predic- 
tions for star-forming galaxies and UHECRs are highly 
model dependent. 

Other astrophysical sources ma y em it in the high latitude 
7-ray sky: i) radio-quiet AGN [33l . |33 | , a nd Fanaroff and 
Riley radio galaxies of type I and II |35l - [37| whose con- 
tribution is strongly model dependent and likely bound 
to few percent of the EGB; ii) 7-ray bursts (GRBs), es- 
timated less than 1% of the diffuse extragalactic 7-ray 
background [38|; iii) star-burst and luminous infrared 
galaxies. The relevant flux may cover a signiflcant frac- 
tion of the EGB (< 20%) [H, but the model dependence 
is such to prevent flrm statements on the relevance of 
this extragalactic source; iv) nearby clusters of galaxies. 



which could yield about 1% - 10% of the EGRET EGB 
[40I - I42I : v) gravitational induced shock waves, produced 
during cluster mergers and large-scale structure forma- 
tion, whose fluxes are quite model dependent and may 
reach few percent [i^, |4J] . All these 7-ray sources have 
been shown to contribute to less than 1% of the Fermi- 
LAT EGB or to be too highly model dependent. In the 
latter case, a very high uncertainty band would be as- 
sociated with the 7-ray source, whose lower limit likely 
gives a negligible contribution to the Fermi-LAT EGB. 
They will therefore be neglected in rest of our paper. 
In the following, we describe few classes of 7-ray emit- 
ters whose unresolved flux is firmly estimated in a non- 
negligible Fermi-LAT EGB percentage. In a conservative 
scenario (Model I), we will subtract AGN and MSPs to 
the Fermi-LAT EGB as derived in Ref. 0. A more 
relaxed model (Model II) will be drawn by the further 
subtraction of a minimal flux from star-forming galaxies 
and CRs at the highest energies. 



1. BL Lacs and FSRQs 

Blazars constitute the class of 7-ray emitters with the 
largest number of identified members. Therefore, unre- 
solved blazars are expected to have a sizable contribution 
to the EGB, [i^. The largest uncertainties in determin- 
ing the blazars contribution are their unknown spectral 
energy distribution and luminosity function (s^ . |46| . In 
addition to phenomenological predictions, an analysis of 
the observed source count distribution through Monte 
Carlo simulations has been performed in Ref. [2J|. The 
reliability of the algorithm relies onto a good agreement 
with the real data, from the comparison of reconstructed 
7-ray fluxes and spectral properties of the sources. The 
energy spectrum is well described by a power- law for both 
FSRQs (softer) and BL Lacs one (harder), being the in- 
tersection between the two fluxes at about 400 MeV. Fol- 
lowing our conservative approach - which is meant to con- 
sider the minimum unavoidable contribution to the EGB 
from unresolved astrophysical sources - we will adopt 
blazar contributions from the curves delimiting the lower 
uncertainty bands displayed in Fig. 20 of Ref. ^29] . The 
ensuing flux is displayed in our Fig. [T] as dotted (dot- 
dashed) line for BL Lacs (FSRQs) contribution. 



2. Pulsars and MSPs 

As a result of their short periods, typical MSPs may 
be brighter in the 7-rays and much older than ordinary 
pulsars [47|] . The ages of MSPs generally exceed the os- 
cillation time across the galactic disk by a large factor 
so that MSPs are expected to be more prevalent at high 
latitudes. On the contrary young, energetic ordinary pul- 
sars are more concentrated close to the Galactic plane, 
where they were born. In the flrst year of Fermi-LAT ob- 
servations 28], 63 pulsars have been identifled. Among 
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them: (i) 16 pulsars at |6| > 10°, of which 11 are MSPs; 
(ii) 5 MSPs at |6| > 40° and (iii) 1 MSP at |6| > 60°. 
We estimate a minimal but not negligible contribution 
of the unresolved MSPs population to the 7-ray flux at 
high latitudes. We adopt an empirical prescription out- 
lined in Ref. [i^, which is based on the spectra of the 
eight MSPs detected by Fermi in the first 9 months 'S] of 
operation. The differential energy spectra of the Fermi- 
detected MSPs are well described by a truncated power 
law: 



dN 
dE 



^ KE 



(1) 



1.9 



r and Ecut are assumed to be (F) = 1.5 and {Ecut) 
GeV, while K has been obtained for \h\ > 40°. 
In order to evaluate Eq. ([1]) for different observational re- 
gions - namely changing the normalization K - we follow 
the prescriptions given in Ref. (47| . Assuming a disk- 
like latitude profile, the ratio of the average intensities at 
different latitudes is given by: 



lMSp{\b\ > h) ^ ln[(sin|biri] 
lMSp{\b\>b2) ln[(sin|b2|)-i]^ 



(2) 



where lMSp{\b\ > bi), i = 1,2, is the average MSP in- 
tensity over a solid angle ft = 4tt{1 — sin 1 6^1) defined by 
the integration from the minimal latitude hi up to 90°, 
written as: 
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n 



1 



s, 



th 



l-<5 



■NiS > Sth). 
(3) 

Smin refers to the assumed Euclidean logN-logS flux 
distribution of the galactic MSP population, which is 
parametrized by a power-law with spectral index 5 — 2.5 
for S > Smin- According to 0, we set Smin = W^^° ph 
s^^ cm^^. A^(> Sth) is the number of resolved sources 
above a given flux threshold Sth- We update the esti- 
mation for I MSP in Ref. [13] with the more recent ob- 
servations for |5| > 10° reported in Ref. [i^, where 8 
MSPs have been found above Sth — 2 ■ 10~® ph s~^ cm"-^ 
(lowest detected MSP flux). We find: 

hisp{\b\ > 10°, E> 100 MeV) = 6.54-10"'^cm"2g-ig^-i 

(4) 

K is then derived from: 



I MSP 



dN 

'dE' 



dE 



(5) 



where ^ refers to Eq. ([T]). Cross-checking the average 
MSP intensities obtained with the prescription outlined 
above, and the results in Ref. 48| for |6| > 40°, we find 
a relative difference of about 30%, due to the theoretical 
uncertainties on the assumed logN-logS and the latitude 
profile. We consider such a discrepancy as an empirical 
theoretical uncertainty on the determination of K , and 
fix the unresolved MSPs contribution subtracting a 30% 
uncertainty from the estimated average intensity. The 
MSP contribution is shown in Fig. [T] as a double dot- 
dashed hue. 



3. Star-forming Galaxies 

Unresolved normal star-forming galaxies are expected 
to give a guaranteed contribution to the high latitude 
isotropic diffuse 7-ray background. Fermi-LAT has iden- 
tified the source of the diffuse emission from our Galaxy 
due to the collisions of CRs with interstellar gas, lead- 
ing to 7-rays from 7r° decay in fiight. This observation 
provides a ground to estimate the 7-ray luminosity of 
star-forming galaxies, by scaling the CR flux with the 
massive star formation rate and fixing the amount of the 
gas in the external galaxy. Theoretical predictions are 
greatly affected by uncertainties in the determination of 
the star formation rate of the galaxies and their gas con- 
tent dll, m, 113] . Given the uncertainty surrounding key 
elements of the determination of this contribution, in our 
strictly conservative approach we do not take into ac- 
count this component. In a more relaxed perspective, 
we consider the lowest predicted contribution from star- 
forming galaxies [3l|. It is derived assuming an increase 
in the number of star forming galaxies with the redshift. 
The adopted emission corresponds to the long dashed 
curve in Fig. [2] 



4. UHECRs 

UHECRs accelerated in astrophysical objects produce 
secondary electromagnetic cascades during their prop- 
agation in the cosmic microwave and infrared back- 
grounds. Ref. [s^l shows that if the primary CRs are 
dominated by protons, such cascades can contribute be- 
tween 1% and 50% of the GeV-TeV diffuse photon fiux 
measured by the EGRET experiment. In Ref. [5l|, the 
EGB spectrum from UHECRs (normalized to the HiRes 
data) has been obtained through a Monte Carlo simula- 
tion of the cascade development and compared with the 
measurement of the EGB by Fermi-LAT. In our more re- 
laxed, whether conservative scenario, we will subtract the 
ankle model contribution to the Fermi-LAT EGB |5l| . 
which we show in Fig. [2l This 7-ray component has the 
peculiar behaviour to slightly increase with increasing en- 
ergy, and at 100 GeV may account 8% of the Fermi-LAT 
measured EGB. 



A. Models for the EGB 

As a result of the previous analysis, we now proceed by 
subtracting from the Fermi-LAT EGB [Tj additional con- 
tributions from unresolved sources at latitudes \h\ > 10°. 
The contributions to the EGB that we will remove from 
the Fermi-LAT spectrum are minimal. In fact, as exten- 
sively explained in the previous sections, the predictions 
that we will take into account for the relevant unresolved 
sources are the lowest ones according to the literature. In 
addition, for MSPs we have lowered existing calculations 
by updating them to the Fermi-LAT observations. 



4 



LU 
T3 









- 


Fermi EGB 






Model 1 EGB 






Sum ol contributions 






FSRQs 






BL Lacs 






MSPs(b>10) 






































t 



10"' 1 10 10^ 

E (GeV) 



FIG. 1: 7-ray spectrum for |6| > 10° latitudes. Fermi-LAT 
data points are displayed along with their power-law fit (solid 
black curve) 7]. The dotted (blue), dot-dashed (green) and 
double dot-dashed (purple) curves correspond to BL Lacs, FS- 
RQs and MSPs contribution, respectively. The dashed (red) 
curve is the sum of the previous three fluxes. The solid (lower, 
green) curve is derived by subtracting the three contributions 
to the Fermi-LAT resuh (Model I). 
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resolved MSPs flux obtained according to the prescrip- 
tion Sect. Ill 21 The results are shown in Fig. [U where 
the Fermi-LAT EGB data are shown along with our 
power-law fit. The contributions from BL Lacs, FSRQs 
and MSPs are identified by dotted, dot-dashed and dou- 
ble dotted-dashed curves, respectively. The fluxes from 
the blazar populations follow power-laws, with softer 
(harder) spectrum for the FSRQs (BL Lacs). The cross- 
ing point for the two curves is around 300 MeV: above 
this energy BL Lacs flux dominates over the FSRQs one. 
The 7-rays from unresolved MSPs show a peculiar spec- 
trum peaked at about 1 GeV and dominate over the 
blazar spectra from 300 MeV up to 3-4 GeV. The sum 
of the three contributions reflects the MSPs flux shape 
with a mild bump. At about 100 MeV the three sources 
explain 10% of the Fermi-LAT EGB and 30% above 1 
GeV. The residual flux Model I, obtained by subtracting 
the sum of the three contributions (dashed curve) to the 
Fermi-LAT best flt flux is identifled by the lower solid 
curve. It is not a net power law due to the dip in the 
GeV region introduced by the MSPs flux. 
Fig. [5] refers to the scenario where the additional contri- 
butions from star-forming galaxies (long dashed line) and 
UHECRs (solid points) as outlined in Sect. Ill 31 and Sect. 
Ill 41 add to explaining the Fermi-LAT EGB. These two 
further contributions add with the previous ones (blazars 
and MSPs of Model I) and the total sum is displayed by 
the dashed red line. The solid (green) curve derives from 
the subtraction of all these contributions from the Fermi- 
LAT EGB (solid black line fltting the data points) and 
is labelled Model II hereafter. Notably, the contribu- 
tion from star-forming galaxies turns out to be relevant 
for < 1 GeV, whereas the 7-rays from UHECRs give 
non-negligible fluxes only at the high-end of the energy 
spectrum. We notice that at 100 MeV Model II explains 
about 70% of the Fermi-LAT EGB, while above 1-2 GeV 
they count about 50% of the total. To consider additional 
astrophysical components to the EGB further decreases 
the residual flux (lower solid line) with respect to the 
Fermi-LAT EGB (upper solid line) and shrinks the room 
left to potential exotic sources, like DM annihilations. 



FIG. 2: 7-ray spectrum for |6| > 10° latitudes. Fermi-LAT 
data points are displayed along with their power-law fit (solid 
black curve) 01 • Dots and long dashed-curve (light blue) 
correspond to the UHECRs and star-forming galaxies 7-ray 
fluxes, respectively. The short-dashed (red) curve corresponds 
to the sum of BL Lacs, FSRQs and MSPs contribution (see 
Fig- HI) the short-dashed (blue) to the sum of the previ- 
ous components with the star-forming galaxies and UHECRs 
ones. The solid (lower, green) curve is derived by subtracting 
all the contributions to the Fermi-LAT result (Model II). 



III. UPPER BOUNDS ON DM ANNIHILATION 
CROSS SECTION 

In this Section we derive conservative upper limits on 
the WIMP annihilation cross section. We make the hy- 
pothesis that the residual fluxes we have derived in Sect. 
Ill Al are entirely provided by the 7-rays produced by ther- 
malized WIMP DM in the halo of the Milky Way. 



A. 7-rays from DM annihilation 

In what we label Model I, we subtract from the Fermi- 
LAT EGB [7] the unresolved contributions for both BL The flux of 7-rays {E^ , ij^) originated from WIMP 
Lacs and FSRQs as outlined in Sect. Ill II and the un- pair annihilation in the galactic halo (52l - [53 | and coming 
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Halo profile 


Isothermal 
a — 3.5 kpc 


NFW [58] 
a = 25 kpc 
rc = 0.01 pc 


Einasto [59] 
a = 0.142 
r_2 = 26.4 pc 
p-2 = 0.05 GeV cm"^ 


|6| > 10° 


2.389 


2.400 


2.833 


10° < l&l < 20° 


4.020 


4.166 


5.752 


|6| > 60° 


1.226 


1.283 


1.232 



TABLE I: Values for Iaq in units of GeV^cm-'^kpc. For all 
these profiles pi — 0.4 GeV cm~"^ , Rsun = 8.2 kpc. 



from the angular direction tp is given by: 



(av) dN^ 1^ 



1 



47r ml dE^ 2' 



(6) 



where {av) is the annihilation cross section times the rel- 
ative velocity mediated over the galactic velocity distri- 
bution function, and dN^/dEj is the energy spectrum of 
7-rays originated from a single DM pair annihilation. In 
particular, we may identify WIMP candidates with neu- 
tralinos in the Minimal Supersymmetric Standard Model 
(see Ref. [HJl and refs. therein). 

The photon spectrum in the continuum originates from 
the production of fermions, gauge bosons, Higgs bosons, 
and gluons from the annihilation of WIMP pairs. The 
spectra dN^/dE^ from DM final states into 66, and 
T+r~ have been taken from Refs. [5^[5^. The extrap- 
olation down to TD.^=\Q GeV seems guaranteed within 
10% of uncertainty for all the annihilation channels [s^l 
(a more careful derivation being beyond the scope of the 
paper). 

The quantity I{ip) is the integral performed along the 
1.0. s. of the squared DM density distribution: 



9'ir{X,i;))dX. 



(7) 



with ■(/; being the angle between the l.o.s. and the di- 
rection pointing toward the galactic center (GC) and 
defined in function of the galactic coordinates so that 
cosV' = cos 6 cos L When comparing with experimental 
data, Eq. ([7]) must be averaged over the telescope ob- 
serving solid angle. Aft: 



1 

An 



I{'ip{b,l))dn. 



(8) 



An 



The integral of the squared DM density over the line-of- 
sight depends from the choice on p{r). When including 
the galactic center in the integration (Eq. (0), different 
DM distributions may lead to very different results for 
/(V'). However, since our analysis is applied to high lat- 
itude regions, the various descriptions for p{r) point to 
very similar values for /(i/'). We neglect any dumpiness 
effects and assume a smooth distribution of DM in the 
galactic halo. The results for /ao for different DM den- 
sity distributions and observational regions are reported 



in Table m All the DM profiles provide very similar results 
for latitudes well above the galactic plane. Hereafter, the 
results will be provided for the cored isothermal density 
profile. 



B. Results on annihilation cross section 
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FIG. 3: Upper bounds on (av) from 7-ray in the high latitude 
galactic halo, as a function of the DM mass. From top to bot- 
tom, solid lines refer to 90% C.L. limits from the comparison 
with Fermi-LAT EGB (black lines), Model I (red lines), Model 
II (blue lines) (see text for details). Dotted, solid and dashed 
lines correspond DM annihilation into p'^ fi~ , bb, t^t~ , re- 
spectively. 

In this Section we derive upper bounds at 90% C.L. 
on the WIMP annihilation cross section from the 7-ray 
Fermi-LAT EGB and the EGB residual fluxes identified 
as Model I and II in Sect. UTJ\ 

For the Fermi-LAT EGB the upper bounds at 90% 
C.L. on {(Tv) are obtained by requiring that the DM sig- 
nal calculated according to Eq. IH] does not exceed the 
measured flux plus 1.28(7 (one-sided upper limit on the 
{av) parameter). The corresponding constraints are plot- 
ted as black lines in Fig. [31 From the same data we have 
then subtracted the unresolved blazars and MSPs min- 
imal contribution, as described at length in Sect. Ill A[ 
and derived the upper bounds on {av) corresponding to 
Model I (red lines). Similarly, upper bounds for the 
Model II EGB are obtained from the further subtrac- 
tion of the minimal flux from star forming galaxies and 
UHECRs (blue lines). In Fig. [3] we display the conserva- 
tive upper bounds on the thermal annihilation DM cross 
section at 90% C.L., derived within the previous assump- 
tions. From top to bottom, each bunch of lines refer to 
the limits on (av) , arising from the comparison with the 
Fermi-LAT 90% C.L. EGB, Model I and Model II. Dot- 
ted, solid and dashed lines correspond DM annihilation 
into fj,'^ fi~ , bb, T~^T~, respectively. Given the scaling of 
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the DM flux oc Ti"^, constraints on {av) increase with 
the mass and span about two orders of magnitude in the 
considered mass interval. It is evident from Fig. |3]that 
the subtraction of the minimal amount of 7-rays from 
unresolved sources lowers the limits on {av) by at least 
50%. The Fermi-LAT data for the EGB are available also 
for latitudes 10° < |6| < 20° and |6| > 60° R. The flux in 
Eq. ([B]) changes for the mere normalization factors given 
in Table m However, given the intensity of the measured 
fluxes our upper limits do not change if derived for the 
other high latitude regions. 

Given the theoretical uncertainties affecting the DM 
content and the astrophysical backgrounds, the results 
in Fig. |3] are of the same order of magnitude or lower 
than the bounds on {av) from cosmological DM 18], 
from the galactic center [26.] . or from inverse Compton 
processes evaluated from 7-rays in different portions of 
the sky [23] ■ Very recently (during the review process 
of the present paper) the Fermi-LAT collaboration per- 
formed a combined analysis on ten Milky Way satellite 
galaxies [l^j, corroborated by the analysis in Ref. [? ]. 
The absence of DM signals from these objects leads to 
upper limits on {av) which are close to 10~^^ for masses 
about 10 GeV and lO-^* for m^=l TeV. These bounds 
are close to the ones established in the present work for 
the high mass side, and stronger for the low mass range. 
The two results, given unavoidable modelling in the ex- 
traction of the upper bounds, strengthen each other in 
disfavouring a DM candidate with an annihilation cross 
section much higher than the electroweak reference value 
3 • 10^^^ cm^/s for very low WIMP masses. We make 
notice that the EGB spectra we have obtained in Model 
I and II could be further reduced, whether by the sub- 
traction of additional components or by increasing the 
predictions for each contribution, set at the minimum in 
the present work. A smaller 7-ray flux at high latitudes 
could therefore be as powerful as the measurements from 
the dwarf spheroidal galaxies. 

We emphasize that our limits arc almost model indepen- 
dent: little dependence on the DM distribution, being at 
high latitudes, and mild differences due to final states. 
Our limits are conservative: it is very unlikely that a 
higher {av) be compatible with Fermi-LAT EGB. Sim- 
ilarly, our upper limits could be lowered only with as- 
sumptions on non-homogeneous DM distributions or, of 
course, comparing to a smaller EGB residual. 



C. Bounds on the Sommerfeld enhancement for 

{(jv) 

Recent claims on the excess of CR positrons [l| have 
stimulated the interpretation of data in terms of anni- 
hilating DM with fairly large annihilation cross sections 
of the order of 10"^^ - lO"^^ 

cir? /s. These numbers 
are at least three orders of magnitude larger than the 
value indicated by observations of the DM abundance 
due to thermal production. One way to boost the an- 



nihilation cross section is through the Sommerfeld effect 
[H, [60l - l6^ , generically due to an attractive force acting 
between two particles, i.e. a Yukawa or a gauge interac- 
tion. In the case of DM particles, the main effect of such 
an attractive force would be to enhance {av) by a fac- 
tor proportional to 1//9 — c/v, where v is the velocity of 
the DM particle [l/v enhancement). The net result on 
the annihilation cross section writes as {av) = S {av)o, 
where S sizes the Sommerfeld enhancement of the anni- 
hilation amplitude. We have evaluated the Sommerfeld 
enhancement S using the approximation of the Yukawa 
potential by the Hulthen potential, for which an analytic 
solution is possible [H, 111] (and checked that the solu- 
tion coincides with the numerical one). The Sommerfeld 
enhancement factor behaves as 1 /v and for very small ve- 
locities it saturates to constant values. Given a the cou- 
pling constant and r7i0 the mass of the new force carrier, 
if the quantity /f^x ' close to the values that make 
the Yukawa potential have zero-energy bound states, the 
enhancement is much larger; indeed, the enhanced cross 
section shows resonances at rriy^ = (n = 1,2,3...), 

which grow as 1/w^, up to the point where they get cut 
off by finite width effects. 

In Fig. 2] we show the Sommerfeld enhanced cross sec- 
tions for a = B — 10^^ and a force carrier of mass 
= 1 GeV (upper curve) and = 90 GeV (lower 
curve). We over- impose the upper bounds obtained in 
the previous Section from the residual EGB Model I and 
Model II and already displayed in Fig. [3l Our results 
show that a Sommerfeld enhancement due to a force car- 
rier of m0 < 1 GeV (a = -^) is strongly excluded by 
Model I and II for the Fermi-LAT EGB data. For a 
massive force carrier (90 GeV) only the resonant peaks 
above the TeV mass are excluded. The result holds for 
/3 = 10~* up to /? = 10~^. Comparable constraints have 
been obtained in [2l|, |63] through the analysis of pertur- 
bations to the CMB angular power spectrum 

Therefore, high latitude 7-ray observations interpreted 
as due to DM annihilation in the Milky Way halo bound 
the Sommerfeld enhancement of the annihilation cross 
section to a factor of 3-10-50-200 for to^ = 10-100-1000- 
5000 GeV, respectively. In case a Yukawa-like potential 
describes this non-relativistic quantum effect, a force car- 
rier heavier than 1 GeV is definitely required. 

D. Bounds from the high-redshift protohalos 

A possible way to boost the annihilation rate is to mod- 
ify the particle theory and make the ansatz that the anni- 
hilation cross section depends on the inverse of the veloc- 
ity [6^. A boosted production of 7-rays in models with 
(crw)oc l/v has been proposed for the first bound objects 
formed in the early phases of the universe [gl, [t^] . After 
the matter-radiation equality is reached, DM perturba- 
tions start growing via gravitational instability and form 
the first bound protohalos at a redshift of about 140. The 
birth of these protohalos depends on the properties of the 
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FIG. 4: Sommerfeld enhancement of the annihilation cross 
section as a function of the DM mass, for a = Solid 

' 47r 

curves are for ^ = 10~^, dotted ones for 13 = 10"^. The 
upper (lower) resonant curve is obtained for a force carrier of 
mass = 1 GeV (90 GeV). The upper (lower) dotted, solid 
and dashed curves correspond to the upper bounds for EGB 
Model I (Model II) derived from WIMPs annihilating in the 
high latitude galactic halo in bb, t^t~ , respectively 

(see Fig. El) 



predicted by: 

p, = 5.28.10^ (^) ' GeVcm-^ 

(10) 

where the cosmological boost factor B, normalized to 
2.6 (i?2.6 — B/2.6) takes into account that the DM is 
distributed according to a Navarro-Frenk- White (NFW) 
density profile with the lowest concentration parameter 
[70| . Eq. (jlOp can be compared with the experimental 
photon density inferred for the Fermi-LAT EGB 7] and 
for our two EGB models derived in Sect. Ill A| which is 
obtained by integrating the photon flux on the Fermi- 
LAT energy range (100 MeV - 100 GeV). We obtain: 

~ 6.62-10-^^ f j GeV cm'^ (Fermi - LAT) 

(11) 

/ c; N -0.41 

p^ ~ 5.65 • 10-^^ f j GcV cm-3 (Model I) 

(12) 



DM particles, since they are responsible of the primor- 
dial inhomogeneities. The complete decoupling of DM 
particles from the thermal bath happens later with re- 
spect to the freeze-out temperature Tf, at a temperature 
of kinetic decoupling Tkd, because scattering events with 
SM particles keep the WIMPs close to thermal equilib- 
rium. For T < Tkd, free-streaming and acoustic oscil- 
lations compete to damp the power spectrum of matter 
density fluctuations, which sets in turn the mass Afc of 
primordial DM structures [7ll-l73|. For very small values 
of Mc, the details of the QCD phase transition could 
further slightly damp the actual cutoff mass [75|. In 



general, however, both the thermal and kinematic de- 
couplings from cosmic plasma are heavily linked to the 
WIMP nature and interactions with SM particles. The 
velocity dispersion of the first protohalos that collapse 
at redshift zc is estimated to be very small {/3 ^ 10"^) 
[70| . Therefore, models for (av) depending on the inverse 
of V predict a boosted flux of DM annihilation products. 
The photons arising from WIMP annihilations in very 
early halos can freely propagate with their energy red- 
shifting and reach the Earth in the range ~ keV - TeV, 
while photons emitted out of this transparency window 
are absorbed by the intergalactic medium. 

The 1 /v enhancement of the annihilation cross section 
may be simply parameterized by writing 70|]: 



{(Tv) = {av)o— cm'^/s. 

V 



(9) 



The energy density in photons today from WIMP an- 
nihilation in the primordial halos can be theoretically 



4.5 • 10 



-16 



E~, 



-0.46 



GeV cm 



GeV, 

(Model II, E.^ > 8 GeV) 



(13) 



We constrain (crw)o by comparison of the theoretical ex- 
pression [TU] with the experimental 7-ray density. The re- 
sults are displayed in Fig. [5] as a function of the WIMP 
mass. The three central lines bound the {av)o (Eq. (|9])) 
parameter from Fermi-LAT photon density given in Eq. 
(ITT]) . Eq. ([T^ and Eq. respectively, from top to 

bottom, when Mc = M^. The upper (dotted) and the 
lower (dashed) bounds are derived for Model II when 
Mc = lO^Af© (upper) and IO'^Mq (lower). The bounds 
on (ctw)o are strong: for WIMP masses below 100 GeV 
it is forced to be < lO^'^'^ cm^/s. Upper bounds grow 
to < 10^32 cm^/s for ~ 1 TeV and sets to < lO'^^ 
cm^/s at 10 TeV. We make notice that they are more 
stringent than limits obtained from primordial light ele- 
ments abundance and CMB anisotropies [tB] and signif- 
icantly improve the bounds of Ref . ^^ . 
The 1/v behaviour of the (av) may be identified with 
the Sommerfeld effect for velocities /? » (m^/m^)^/^ 
j6ll |. For lower velocities, as are the ones typical for pro- 
tohalos, the series of resonances appears (see Sect. IIII Cp 
and the Sommerfeld enhancement S behaves as 1 /w^ close 
to the peaks. In this case, the upper bounds on the 
annihilation cross section may be obtained by rescaling 
{av) — {av)QS with a factor 1//3 • m,p/m^. From Eqs. [10] 
- [TT] it is straightforward to notice that the bounds on 
a Sommerfeld enhanced {av) derived from a overproduc- 
tion of 7-rays in protohalos, are much weaker than the 
ones imposed by annihilation in the high-latitude galactic 
halo. 
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(GeV) 

FIG. 5: Bounds on ((t?;)o from Eq. ((9]), as a function of 
the DM mass. The central three bounds are obtained for 
Mc = M®, and from Eqs. HI} (black hue), ((12]) (red line) and 
(|14|l (blue line) respectively, from top to bottom. The upper 
(lower) purple lines are derived from Eq. (|14|) for Model II 
EGB and Mc = 10^M©(10"^M®). 



IV. CONCLUSIONS 

The 7-ray EGB measured by Fermi-LAT [Tij likely 
includes contributions from galactic and extragalactic 
unresolved sources. We have explored possible non- 
negligible diffuse contributions from unresolved blazars, 
MSPs, star-forming galaxies and UHECRs. Lead by 
a conservative attitude, we have considered the mini- 
mal contribution for all sources, and neglected those ob- 
jects whose high latitude flux is not excluded to be less 
than 1% of Fermi-LAT EGB. Two residual EGB fluxes 
have been derived by subtraction of the additional fluxes 
from the Fermi-LAT EGB: Model I is obtained after the 
subtraction of unresolved BL Lacs, FSRQs and galactic 
MSPs, while Model II is the residual flux after the fur- 
ther subtraction of star-forming galaxies and UHECRs. 
From our new residual EGB fluxes, we have set upper 
limits on the DM annihilation cross section into 7-rays. 
A conservative upper bound on {av) is derived by as- 
suming that the Model I and II EGB are entirely due 
to WIMPs pair- annihilating in the halo of our Galaxy. 
Values for {av) > 10^^^ cm'^/s are strongly excluded for 
~ 10 GeV,'~^hile for ~ 100 GeV (1 TeV) the 
annihilation rate is bounded to 3 ■ lO^'^^ cm^/s (10^^'* 
cm'^/s). This results holds for DM annihilating into bb. 
Stronger limits below m^— 1 TeV are derived for annihi- 
lation into the Icptonic r annihilating channel, while for 



the /i channel the limits are close to the bb below m^= 
100 GeV, and weaker above this mass. Annihilation into 
leptons is therefore excluded at a level which strongly 
disfavours the interpretation of cosmic positron fraction 
data in terms of leptophilic DM with small cosmologi- 
cal boost factors. The latter boost factors are in turns 
strongly limited by antiproton data [s^] . 
The bounds on {av) have been interpreted in terms of 
Sommerfeld enhancement of the annihilation cross sec- 
tion. A Sommerfeld enhancement due to a force car- 
rier of m0 < 1 GeV {a — j^) is strongly excluded by 
Model I and II for the Fermi-LAT EGB data. For a 
massive force carrier (90 GeV) only the resonant peaks 
above the TeV mass are excluded. High latitude 7-ray 
observations interpreted as due to DM annihilation in the 
Milky Way halo bound the Sommerfeld enhancement of 
the annihilation cross section to a factor of 3-10-50-200 
for m^=lG-lGG-lGGG-5000 GeV, respectively, and in case 
an annihilation into light quarks occurs. For < 6-700 
GeV these limits are reduced by a factor of few for the 
pure T'^T~ annihilation channel. In case a Yukawa-like 
potential describes this non-relativistic quantum effect, a 
force carrier heavier than 1 GeV is definitely required. 
Finally, we have explored the possibility that the residual 
7-ray EGB is entirely due to cosmological annihilation 
of DM in protohalos at high redshift. Within the hy- 
pothesis that {av) is inversely proportional to the WIMP 
velocity, very severe limits are derived for the velocity- 
independent part of the annihilation cross section, de- 
pending on the protohalo mass. 
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